Pattern formation and phase transitions in systems with non-monotonic interaction 
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We analyzed pattern formation and identified different phases in a system of particles interacting 
through a non-monotonic short-range repulsive (r < re) and long-range attractive (r > re) potential, 
using molecular-dynamics simulations. Depending on parameters, the interaction potential models 
the interparticle interaction in various physical systems ranging from atoms, molecules and colloids 
to neutron stars and vortices in low k type-II superconductors and in recently discovered "type-1.5" 
superconductors. We constructed the phase diagram "critical radius re — density n" and proposed a 
new approach to characterization of the phases. Namely, we elaborated a set of quantitative criteria 
in order to identify the different phases using the Radial Distribution Function (RDF), the local 
density function, and the occupation factor. 



I. INTRODUCTION 

Pattern formation in many systems is governed by 
competing interaction^. Examples of such systems are: 
the pasta phase in neutron stars^, ferrofluids'^"^ , Lang- 
muir monolayers''^, magnetic garnet thin films, type-I 
superconductors^, colloids and gels^""^'*, etc. In general, 
there is a strong correlation between the pattern forma- 
tion and the inter-particle interaction. Attraction favors 
aggregation, while repulsion favors low local densities. 
This competition between repulsive and attractive inter- 
action leads to very rich phases, such as stripes, clusters, 
bubbles, etc.^^. Those different phases were observed in 
many diverse systems. In neutron stars, the competition 
between the short-range nuclear attraction and the long- 
range Coulomb repulsion leads to complex pasta phases^. 
In ferrofiuid systems, rich phases due to the competi- 
tion between dipolar forces and short-range forces op- 
posing density variations were found experimentally"^ and 
theoretically'*. 

Pattern formation was extensively studied in colloidal 
systems where the colloid-colloid interaction is charac- 
terized by the competition of the hardcore excluded 
volume interaction, on the one hand, and the polar- 
ization of the particles, on the other hand. The in- 
teraction potential can be further controlled by adding 
other components. As a result, rich configurations, 
such as clusters, repulsive or attractive glassy states, 
gels, were found numerically-'^^"^"'^ , analytically^^'^"^ and 
experimentally*^'*"^ in colloidal systems with short range 
attractive interaction. The properties of isolated clusters 
formed by the short-range attractive and long-range re- 
pulsive interaction were studied^**. By controlling the 
interaction, the growth of the cluster was shown to 
change from nearly spherical to one-dimensional. The 
one-dimensional growth of the clusters facilitated the col- 
lective packing into columnar or lamellar phases^^. The 
columnar and lamellar phases in three dimension were 
also studied by molecular dynamics (MD) simulations^^. 

In superconductors, the vortex-vortex interaction is 
usually considered to be either repulsive (in type-II su- 
perconductors where the Ginzburg-Landau parameter k. 



i.e., the ratio of the magnetic field penetration depth A 
to the coherence length ^, k > l/-\/2) or attractive (in 
type-I superconductors, where k < l/-\/2 and vortices are 
unstable) while vortices do not interact with each other 
at the so-called "dual point" when k = l/-\/2^^- How- 
ever, a deeper analysis of the inter-vortex interaction in 
type-II superconductors near the dual point revealed an 
attractive tail^^'^*. This repulsive-attractive inter- vortex 
interaction was used for the explanation of unusual pat- 
terns in the intermediate state in 1ow-k superconductors 
(e.g., Nb): islands of Meissner phase surrounded by vor- 
tex phase or vice versa, i.e., vortex clusters surrounded 
by Meissner phase^^'^*. 

Recent discovery of "type-1.5" superconductors^^ in- 
duced a new wave of interest (see, e.g., Ref."^"'"^*) to 
systems with non-monotonic interactions, due to the 
fact that the observed vortex patterns in those super- 
conductors revealed a clear signature of the repulsive- 
attractive inter-vortex interaction. In particular, several 
properties (e.g., vortex lattice with voids, the nearest- 
neighbor distribution) of the observed vortex patterns 
were explained using a simple model that involved a non- 
monotonic inter-vortex interaction based on a more gen- 
eral approach to multi-ordcr-paramcter condensates'^^. 

In this paper, we consider a model competing inter- 
action potential which is repulsive for short range and 
attractive for long range. This form of the interaction 
potential could be used as a model for different systems 
with non-monotonic intcr-particlc intcration, e.g., atoms 
or molecules (i.e., the Lennard- Jones potential), or vor- 
tices in two-band superconductors, depending on specific 
parameters of the potential. We study pattern formation 
in two-dimensional systems, for different interaction po- 
tential profiles, i.e., we distinguish "soft-core" and "hard- 
core" interactions and analyze the transitions between 
different phases. Based on this, we construct a phase 
diagram for different interaction parameters and particle 
densities. We propose a new approach to characterize the 
different phases: instead of qualitative characterization 
of the phases (e.g., clusters or labyrinths), we introduce 
a number of quantitative criteria to distinguish these. 
In particular, different phases are analyzed in terms of 



the Radial Distribution Function (RDF) and additional 
quantities chracterizing, e.g., the local density of particles 
in clusters. 

The paper is organized as follows. In Sec. II, we de- 
scribe the model. The pattern formation for different in- 
teraction parameters is discussed in Sec. III. In Sec. IV, 
we analyze different phases using the RDF and discuss 
criteria for identification of different patterns. The con- 
clusions are given in Sec. V. 



II. MODEL 

The inter-particle interaction potential is taken to be 
of the following form: 

V,, =Vo[^* Ko{b * n,/X) - Ko{n,/X)^ . (1) 

Here, Kq is a first order Bessel function, r^j =| r^ — r^ | 
is the inter-vortex distance, Vq and A are the units of 
energy and length, respectively. [Note that in case of a 
type-II superconductor appropriate units of length A and 
the energy Vq, correspondingly, are the magnetic field 
penetration depth A and Vq = ^q/Stt'^X^, where $o = 
hc/2e (see, e.g., Ref.^^)]. 

In dimensionless form, the interaction potential (1) 
reads as 



K; -^ = ^*^o(&*r^)-/fo(r:,), 



(2) 



where the dimensionless length is defined as r^ — rij/X. 
Further on, we will omit the primes and use the dimen- 
sionless form of the potential (2). The interaction force 
is then given by 



Fjj = -Wij = (a * Ki{h * Tij) - Ki{r,j))r, 



ij, 



(3) 



where Ki is a first order Bessel function, Yij ~ (r^ — 
Yj)/rij, a and b arc two positive coefficients. 

The interaction potential Eq. (2) is generic: by choos- 
ing suitable parameters, it can be used as a model of 
non-monotonic interactions in different systems, for ex- 
ample, the well-known inter-atomic (molecule) Lennard 
Jones (LJ) potential: fiLjir) = /io[(r/cr)~^^ — (r/tr)"^]. 
In Fig. 1(a), we compare the LJ potential for a particle of 
diameter a — 2.762A and the energy unit ^o = O.IVq with 
the model potential given by Eq. (2) with the parameters: 
a — 1.045x10^, b ~ 5.896. The corresponding interaction 
forces are presented in Fig. 1(b). The comparison shows 
that the potentials and the corresponding forces (i.e., for 
the model potential (1) and the LJ potential) fairly agree 
with each other. On the other hand, one can easily see 
that the interaction potential Eq. (1) is a generalized 
form of the inter-vortex interaction in type-II supercon- 
ductor which, as shown by Kramer^^, can be presented 
in the form: V{r) = di{K)Ko{r/X) - d2iK)Koir/0- 

The first term of Eq. (3) is repulsive while the second 
term describes an attractive interaction force. Indeed, 



for r ^ oo, Ki{x) — > \J tx j {2x)e~'^ . The interaction force 
(3) has a repulsive (attractive) tail when 6 < 1 (6 > 1), 
while for r ^ 0, Kx{x) -^ 1/x and thus F^j -^ [a/b — 
l)/r. Therefore, for short range the interaction force (3) 
is repulsive (attractive) when a > b (a < b), and we only 
consider the case a > b since an attractive interaction for 
short distances would result in a collapse of our system 
of point particles. When a > b and 5 < 1, the interaction 
is always repulsive, and particles form a Wigner crystal 
structure. The most interesting case is realized when a > 
b and 5 > 1, and the interaction has a repulsive core and 
attractive tail. In this case, there exists a critical distance 
Tc, where the inter-particle interaction energy (1) reaches 
a minimum (and the interaction force (3) changes sign). 
By setting the force equal to zero, the coefficient a is 
given by 



a = 



Ki{b*rc)' 



(4) 



The pattern formation is determined by the coefficients 
6, re, and the particle density n. 

We study pattern formation in a system of interacting 
particles using Langevin equations. The dynamics of a 
single particle i obeys the following overdamped equation 
of motion: 



r/v, 



= F, = ^F, 



(5) 



Here, 77 is the viscosity, which is set to unity. F^ is the 
interparticle interaction force defined by Eq. (3) and Fj 
is the stochastic thermal force. 

The thermal stochastic term F^ in Eq. (5) obeys the 
following conditions: 



and 



(^f(t))=0 (6) 



{F^{t)F^{t')) = 2rjkBTS,j6it - t'). (7) 



We consider a two-dimensional (2D) square simulation 
region Lx x Ly in the xy-plane and apply periodic bound- 
ary conditions in the x and y directions. For the interac- 
tion force given by Eq. (3), which exponentially decays 
for large distances, we use a cut-off for r > 8. The length 
of the square cell L^ = Ly = L has been varied from 
60 to 180 to examine the finite-size effects, and we set 
L — 120 to optimize the calculation speed without in- 
fluencing the results. To obtain stable particle patterns, 
we performed simulated annealing simulations (SAS) of 
interacting particles. For this purpose, particles were ini- 
tially randomly distributed inside the simulation region 
at some suitable non-zero temperature (depending on the 
inter-particle interaction). Then temperature was grad- 
ually reduced to zero, and the simulation was continued 
until a stable state was reached. 
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FIG. 1. (Color online) The LJ potential and the model 
potential given by Eq. (1) with the following parameters: 
a = 1.045 X 10^ h = 5.896, a = 2.762A, and ^o = O-l'K) 
(a). The corresponding interaction forces (b). 
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FIG. 2. (Color online) The profile of the inter-particle inter- 
action force versus the distance. The main panel and the inset 
panel show the change in the force profile due to an increase 
of Tc and 6, respectively. 



III. PATTERNS AND PHASES 

In this section we will study pattern formation and 
identify different phases depending on parameters of the 
interaction rj,, 6 and the particle density n. In Fig. 2, we 
illustrate the change of the interaction force profile due to 
the increase of Tc and h. On the one hand, the interaction 
force is very sensitive to Tc, since its increase directly 
leads to the increase of the repulsive part in Eq. (3) (see 
the main panel of Fig. 2) . On the other hand, increase of h 



leads to the increase of both the repulsive and attractive 
parts (see the inset of Fig. 2). However, the repulsive 
interaction increases much faster for r < r^. than the 
attractive interaction for r > re- 

Thus, with increasing b it becomes hard to reduce the 
inter-particle distance below r,.. In other words, the in- 
crease of b results in hardening of the core in the interac- 
tion force (3), i.e., the interaction changes continuously 
from a "soft-core" to a "hard-core" interaction. 



A. Soft-core interaction 

To analyze the soft-core regime, we set the coefficient 
h = 1.1. We define the critical density as ric — 8r^^/\/3 
which is the density of an ideal (hexagonal) Wigner crys- 
tal with the lattice constant a = r^ The density is de- 
fined as n = N/S, where S — 120 x 120 is the area of 
the simulation region and N is the number of particles 
(further on, when analyzing different patterns, we will 
refer to either the number of particles N in the simula- 
tion cell or to the density which is n = A^/14400). For 
the case of n < ric, which is considered here, the Wigner 
crystal is not stable due to the attractive interaction. In 
Fig. 3, we show patterns formed by A^ = 1500 particles 
when Tc increases from 1.9 (a) to 2.9 (f). Note that the 
condition n < ric is always fulfilled for all the values of r^ 
in this range. We found that for r^ < 2.1 particles form 
clusters similar to the formation of the "clump" phase 
found in Ref i^^^** (see Figs. 3(a) and (b)). The main 
difference to the patterns found in Ref.^^"^^ is that the 
relatively softer core in our case is compressed due to the 
attractive interaction, and the clusters acquire a circular 
shape. In addition, the interaction between the clusters 
(decaying exponentially for long distances) becomes neg- 
ligible for inter-cluster distance of the order of few to 
several rg. Therefore, the clusters can be considered as 
non-interacting for low densities (although still they do 
not approach each other), contrary to the situation of 
Ref ^^"^^ when a super-lattice is formed due to the long- 
range cluster-cluster repulsion. When re increases, the 
clusters expand. In particular, for Vc > 2.1, the clusters 
start to elongate which is an indication of the instability 
with respect to the transition to the stripe phase. For 
2.1 < Tc < 2.3, a mixed state with both stripes and 
clusters is observed (see Fig. 3(c)). Further increase of 
Tc gradually destroys the cluster phase and leads to the 
formation of the labyrinth phase (see Figs. 3(d) to (f)). 

In order to investigate the influence of the density, we 
gradually increase the number of particles from 1500 to 
10500 in our simulation cell. First, in Fig. 4, we present 
patterns formed by A^ = 5500 particles for varying re. 
Note that n > n^ for N — 5500. As compared to the 
lower density case shown in Fig. 3, in the cluster phase, 
the additional particles lead to the expansion of the clus- 
ters (see Figs. 4(a) and (b)). A mixture of clusters and 
stripes are formed when the additional particles form 
bridges connecting the clusters (see Fig. 4(c)) which will 



be discussed in detail below. For the labyrinth phase 
(shown in Figs. 3(d) to (f) for N = 1500), the addi- 
tional particles fill the empty regions (voids) , resulting in 
the formation of the triangular lattice with varying local 
density (Fig. 4(d)) and, finally, a regular triangular lat- 
tice (see Figs. 4(e) and (f)). Note that the lattice with 
varying local density (Fig. 4(d)) is not stable in systems 
with pure repulsive interaction such as vortices in type-II 
superconductors. 

We analyzed in detail the intermediate regime (i.e., 
corresponding to the transition from clusters to stripes, 
see Fig. 3(c)) where re « 2.3. The patterns for vary- 
ing density are shown in Fig. 5. Fig. 5(a) displays a 
configuration at low density with N — 1500 when many 
individual clusters are formed. When increasing the den- 
sity (see Fig. 5(b)), the clusters connect with each other 
and form a stripe phase. A further increase of the num- 
ber of particles (Fig. 5(c)) results in the formation of a 
mixed phase of interconnected stripes with voids. A very 
interesting and counter-intuitive evolution is observed 
when the number of particles increases from N — 5500 
to iV = 6500 (see Fig. 5(d)): in contrast to the grad- 
ual transition from a void-rich configuration to a lattice- 
rich configuration when N changes from N = 3500(b) to 
N = 5500(c), in case of iV = 6500 we observe a "reen- 
trant" behavior, i.e., the void-rich phase starts to recover 
which is compensated by an increase of the local density 
in the stripes. However, further increasing density re- 
sults in the expansion of the stripes to the empty regions 
which is accompanied by a decrease of the local density 
in the stripes, as shown in Fig. 5(e). The distribution 
of particles becomes more uniform, with only few small 
voids. Finally, for N = 9500 (Fig. 5(f)), we obtain a de- 
formed triangular lattice characterized by a varying local 
density with only one small void. 

Our calculations show that the obtained phases are 
very sensitive to variations in Tc- Thus, if Tj, slightly 
decreases (e.g., r^ — 2.25), the number of clusters greatly 
increases as compared to the case re = 2.3, for the same 
density of particles. For re — 2.25, we also observe the 
transition from a void-rich configuration to a lattice-rich 
configuration. However, since the decrease of re increases 
the attractive component of the inter-particle interaction 
this occurs at much higher density. For even smaller 
values of r^, i.e., r^ < 2.1, we do not observe the lattice 
phase, even for extremely large number of particles (up 
to iV = 20000). 



B. Phase Diagram 

Based on the above analysis of the patterns and phases, 
we constructed a phase diagram in the plane of r^, and 
the number density n (see Fig. 6). For extremely low 
density (n < 0.03), particles form small clusters which 
are well separated, and the patterns are rather insensi- 
tive to variations of Tj,. However, phases become richer 
when density n increases (i.e., n > 0.1). Low values of 




FIG. 3. (Color online) Different patterns for A*' — 1500 parti- 
cles in a unit cell L x L with L — 120 for varying r^: re = 1.9 
(a), 2.1 (b), 2.3 (c), 2.5 (d), 2.7 (e), 2.9 (f). 



re (i.e., re < 2.1) still favor cluster formation for a broad 
range of densities although the clusters become denser 
for large n. With increasing re, clusters become unsta- 
ble with respect to the formation of "bridges" between 
separate clusters which is a precursor of the formation 
of stripes. Stripes are formed in a rather narrow range 
of Tc when re > 2.1 (see Fig. 6). The stripe phase is 
represented by two sub-phases, I and II, i.e., void-rich 
phase (see Fig. 6(i)) and void-poor (lattice- rich) phase 
(Fig. 6(j)). For larger Tj, and n < n^, particles form 
labyrinth structures. However, when increasing the den- 
sity, additional particles fill the empty regions and finally 
they form a deformed triangular lattice with varying den- 
sity. Depending on r^ and n < rie deformed lattice is 
characterized by a varying local density or by the ap- 
pearance of voids. Correspondingly, we distinguish two 
regions (1 and 2 in Fig. 6). Note that in the vicinity of 
the phase boundaries patterns are always mixtures of the 
two phases (e.g., clusters and stripes) except the phase 
transition from the stripe phase to the deformed trian- 
gular lattice. Near this phase boundary, particles either 
form stripes with high local density or deformed triangu- 
lar lattices with lower local density. The probability of 
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FIG. 4. (Color online) Patterns for A^ — 5500 particles in a 
unit cell L x L with L = 120 for varying re', re = 1.9 (a), 2.1 
(b), 2.3 (c), 2.5 (d), 2.7 (e), 2.9 (f). 



FIG. 5. (Color online) Patterns for r^ — 2.3 and varying 
number of particles N (in the computational unit cell L x L 
with L = 120): N = 1500 (a), 3500 (b), 5500 (c), 6500 (d), 
7500 (e), 9500 (f). 



finding the lattice configuration greatly increases when Vc 
or n increases. Finally, for very large density, when the 
average interparticle distance becomes larger than Tc the 
repulsive interaction prevails resulting in the formation 
of an almost perfect triangular lattice. 

The phases presented in Fig. 6 are found in many real 
systems. For example, the cluster phase is found in such 
systems as colloids^"'^^'^^'^'^ and neutron stars'^. The ob- 
tained stripe patters are very similar to those in Lang- 
muir monolayers*. 

Several of the calculated phases were found in super- 
conductors, e.g., the stripe phase in the mixed state of 
type-I superconductors, clusters of Meissner phase or vor- 
tex clusters in the intermediate state in 1ow-k type-II su- 
perconductors^^'^*, the labyrinth phase (i.e., vortex lat- 
tice with voids) or vortex clusters in type-1.5 supercon- 
ductors^^ . In spite of the variety of the physical systems 
and length scales ranging from nano- and micro-objects 
to cosmic objects, the common feature of all these sys- 
tems is a competing attractive-repulsive interparticle in- 
teraction which allows analyzing them within the same 
approach. 



C. Hard-core interaction 

Let us now analyze the influence of coefficient b on the 
pattern formation. As mentioned above (see Fig. 2), an 
increase in b changes the potential from a soft-core to a 
hard-core. Hardening the repulsive core in the interac- 
tion potential generally leads to a decreasing compress- 
ibility of the inner parts of the patterns where particles 
are closely packed. Correspondingly, the patterns change 
as compared to the soft-core regime. 

As shown in Fig. 7, for 6 = 4 and low density {N = 
1500), all the patterns are clusters of different shape. 
Thus for Tc — 1.15 and r^, — 1.3, the clusters are of circu- 
lar shape similar to those found in the soft-core regime 
although for re — 1.3 some clusters composed of a small 
number of particles have polygon shapes. With increas- 
ing Tc, polygon shaped clusters become more favorable, 
which allows us to identify them as short stripes. For 
even larger r^ the clusters become much larger. They 
represent separate islands of triangular lattices. 

For higher density {N = 6500), the variety of phases 
is much richer. Although for r^ = 1.15 circular-shape 
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FIG. 6. (Color online) The phase diagram in the plane "critical radius re — density n" (n) and representative patterns (a to 
m) for the difTerent phases. For extremely low density (n < 0.03), particles form small well separated clusters. For Vc < 2.1, 
particles form clusters. For larger re, the size of the clusters increases (see the change from (d) to (e) and (b) to (c)). The 
increase of the density results in increasing the size of the clusters (see the change from (d) to (b) and (e) to (c) ). Further 
increase of re leads to elongation of the clusters and the formation of "bridges" between them (a). Thus, the configurations 
change gradually from the cluster phase to the stripe phase (i). When re or n become even larger, stripes interconnect and 
form patterns with voids (j). Thus the stripe phase is divided into two sub-phases, I and II. When re become even larger, 
labyrinth-like configurations are formed for n < Ue (see (g) and (f)). With further increase of re or n, deformed triangular 
lattice with (k) or without voids (1) are formed. The deformation of a triangular lattice is reduced for even larger values of re 
and n (m). 



clusters are still observed (see Fig. 8(a)), which grow 
in size for Vc — 1.3 and slightly deviate in shape from 
circular (Fig. 8(b)), for larger r^ = 1.45 deviations of 
the cluster shape from circular become more pronounced 
(see Fig. 8(c)). With further increasing r^ i.e., r^ = 1.6, 
stripes are formed (see Fig. 8(d)) followed by lattices with 
voids for r^ = 1.75 and r^ ~ 1.9. In the hard core regime, 
deformations of lattices occur via the appearance of voids 
instead of varying local density (in the soft-core regime) . 



IV. ANALYSIS OF THE PATTERNS 
A. Radial Distribution Function 

Although different patterns and phases studied above 
are qualitatively well distinguished, it is highly desirable 
to build a set of solid criteria which would allow us to 
identify them in a quantitative manner. For this purpose. 
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FIG. 7. (Color online) Patterns for iV = 1500 particles in 
a unit cell L x L with L — 120, 6 = 4 and for varying r^: 
Tc = 1.15 (a), 1.3 (b), 1.45 (c), 1.6 (d), 1.75 (e), 1.9 (f). 



FIG. 8. (Color online) Patterns for A^ — 6500 particles in 
a unit cell L x L with L = 120, 6 = 4 and for varying Vc'. 
Tc = 1.15 (a), 1.3 (b), 1.45 (c), 1.6 (d), 1.75 (e), 1.9 (f). 



we analyze here different phases by means of the Radial 
Distribution Function (RDF). The RDF, g{r), describes 
the variation of the atomic (particle) density as a function 
of the distance from one particular atom (particle) . If we 
define an average density as n = N/V (where V is the 
volume or surface area in the 2D case), then the local 
density at distance r is ng(r). The knowledge of RDF 
is important since one can measure g{r) experimentally 
using neutron scattering or x-ray diffraction scattering^^. 
Moreover, macroscopic thermodynamic quantities can be 
calculated using g(r) in thermodynamics'^^. 

In our calculations, we define the RDF gi{r) as follows: 



9i{r) 



AN/Ar 
2'Kr 



(8) 



Here, the lower index indicates that the RDF centers at 
the position of the ith particle; AN is the number of 
particles whose distance to the zth particles is between r 
and r + Ar. The average g{r) is given by 



ffW 



1 ^ 



(9) 



In Fig. 9(a), we plot the function g{r) calculated for the 
low-density cluster configuration shown in Fig. 3(b). The 
function g(r) has two well-pronounced peaks. The first 
peak corresponds to the average distance to the first co- 
ordination sphere (nearest neighbors) while the second 
one is located at the distance approximately twice the 
distance to the first peak, which shows short-range pe- 
riodicity. For r larger than the size of the cluster, and 
smaller than the inter-cluster distance, g(r) < 1, since 
only few particles are located inside this range. As well 
as for the RDF obtained for the pasta phase in neutron 
stars^ , the decreasing tail in our case also shows a strong 
aggregation of the particles in the cluster phases. The 
minimum of g(r) is very close to zero since most of the 
clusters have circular symmetry and they are well sep- 
arated. The position of the minimum of g{r) (marked 
by the gray arrow in Fig. 9(a)) gives an estimate of the 
average diameter of the clusters. 

For the stripe phase (see Fig. 9(b) for the pattern 
shown in Fig. 3(d)), there are also two peaks indicating 
the short-range periodicity, similar to the cluster phase. 
However, unlike in the case of clusters the minimum of 
g(r) is non-zero since the stripes are generally not sep- 
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FIG. 9. (Color online) The average RDF of the patterns 
formed at low density (A'^ — 1500). (a), (b), and (c) cor- 
respond to the patterns shown in Figs. 3(b), (d), and (f), re- 
spectively. The peaks marked by dark arrows show the shot- 
range periodicity. The gray arrow shows the minimum of the 
RDF. 



arated. For the labyrinth phase (see Fig. 9(c) for the 
pattern shown in Fig. 3(f)), only short-range periodic 
ordering exists, and g{r) becomes uniform for larger r. 

In Fig. 10, we plot g{r) for the patterns formed at high 
density (N=5500). For clusters, the increase in the den- 
sity leads to an increasing variation of the cluster size. 
Note, however, that the sizes of individual clusters are 
rather insensitive to moderate variations of the number 
of particles in the clusters, due to the strong compress- 
ibility of the core. As a result, the function g{r) for the 
high density clusters misses the short-range periodicity 
and the second peak disappears (see Fig. 10(a)). Strik- 
ingly, this peak still exist for the void-rich phase (see 
Fig. 10(b)). The decrease of the tail of ^(r) becomes very 
slow, which shows a minor aggregation of the particles. 

For the lattice phase (see Fig. 10(b)), the position of 

the first peak is ri sa a = i/2/(v3n), where a is the 
distance between two neighboring particles in the ideal 
triangular lattice with the density n. The positions of the 
second and third peaks are at r2 ~ Via and at r2 ~ 2a, 
respectively, corresponding to those in a triangular lat- 
tice. A fourth, fifth and even sixth peaks appear, which 
shows that this phase is much more ordered. However, 
due to the variation of the local density, these peaks are 
strongly broadened, and they are actually a combination 
of many neighboring peaks. These peaks become clearer 
for larger Vc or n, which implies that the lattices become 
more regular. 

It is interesting to compare the results of the function 
g{r) for hard-core interaction with those for soft-core in- 
teraction. Although the clusters shown in Figs. 8(a) and 
(b) have shapes very similar to those of the clusters shown 
in Figs. 4(a) and (b), the analysis using the RDF shows 
that most of the clusters shown in Fig. 8(b) have hexag- 
onal ordered cores. The peaks in g(r) appear to be much 
better separated than for, e.g., the deformed triangular 
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FIG. 10. (Color online) The average RDF of the patterns 
formed at high density (N=5500). (a), (b), and (c) correspond 
to the patterns shown in Figs. 4(b), (d), and (f), respectively. 
The arrows have the same meaning as in Fig. 9. 




FIG. 11. (Color online) The average RDF of the patterns 
formed in the hard-core case, (a), (b), (c), (d), and (e) cor- 
respond to the patterns (a), (b), (c), (d), and (e) in Fig. 8, 
respectively. 



lattice formed in the soft-core case (see Fig. 10(c)). For 
the stripe and void-rich phases, the RDF shows much 
clearer triangular lattice ordering. 

Therefore, the analysis using the RDF allowed us to es- 
tablish quantitative criteria for the different phases and 
reveal the differences in the structure of the patterns 
(which look similar, e.g., clusters) in case of soft- and 
hard-core interparticle interaction. 
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FIG. 12. (Color online) The average local density I versus 
the particle density n for b — 1.1 (a). Curves A, B, C, D, 
E, and F correspond to r^ = 1.9, 2.1, 2.3, 2.5, 2.7, and 2.9, 
respectively. For the cluster phase, / > 20, for the stripe 
phase, 10 < 7 < 20, for the labyrinth phase, 6 < 7 < 10. 
7 = 6 for lattice and deformed lattice configurations. The 
local density 7 versus the number of particles A'^ for b — 4 (b). 
Curves A, B, C, D, E, and F correspond to re = 1.15, 1.3, 
1.45, 1.6, 1.75, and 1.9, respectively. 



/ since only the "shell" of those regions of thickness £, is 
considered. In Fig. 12(a), we plot the function I for dif- 
ferent Tc and N for the soft-core interaction with b = 1.1. 
We see that the function / can be used to charaterize 
the different phases. Thus we found that, for clusters, 
/ is always large (/ > 20) since aggregation is strong. 
For stripes, aggregation is smaller, and thus / is smaller 
(10 < I < 17). For labyrinths, the regions free of parti- 
cles are relatively small. Therefore, / ranges between 6 to 
9. Finally, for the lattice phase, / is always 6. Therefore, 
the function / serves as a measure of the aggregation of 
particles and allows us to effectively distinguish the dif- 
ferent phases. 

The function / also provides a tool to analyze the sta- 
bility of patterns with increasing density. This is demon- 
strated in Fig. 12(b), where we plot / for the hard-core 
case when 6 = 4. For Vc = 1.15, / > 20 for all the values 
of the density, and thus the clusters are stable. However, 
the situation changes for re = 1-3: while clusters are 
formed for low densities up to iV « 7500, for larger den- 
sity / decreases below the value / = 20 which means that 
clusters elongate and interconnect giving rise to the onset 
of the stripe phase. For re = 1.45, the particles start to 
form stripes at even lower density (N=3500). For higher 
density, the additional particles fill in the empty regions 
and finally they form the lattice phase. For re > 1.6, 
the island-like clusters formed at low densities are very 
unstable with respect to an increase in density. The ad- 
ditional particles rapidly occupy the empty regions, and 
the patterns change form clusters to stripes and from 
stipes to lattice with increasing density. 

Therefore, the phenomenological description of pattern 
evolution with increasing density we revealed above in 
Sec. III. A has been verified in a quantitative manner us- 
ing the rigid basis of the local density function approach. 



C. Occupation factor 



B. Local density 

Let us define, using the RDF, an "order parameter" to 
characterize the different phases. By definition, the local 
density is: 



2-Krngi{r)dr = N^ — 1. 



(10) 



Here, N^ ~ Tr^^n is the average number of particles 
within the circle centered at the ith particle with ra- 
dius R = ^. Since in an ideal triangular lattice one par- 
ticle has six nearest neighbors, we take N^ = 7, then 
£, = y^N^/nn — ^Jljwn. Thus, for any configuration 
characterized by a small local density fluctuation, the 
average local density / — Iji) — 6. From the defini- 
tion given by Eq. (10), we can see that the presence of a 
large fraction of empty regions can considerably increase 



We demonstrated that the (average) RDF g(r) 
(Eq. (9)) and the local density function / (Eq. (10)) 
allowed us to unambiguously characterize the different 
patterns and phases. At the same time, it is still hard to 
distinguish, using the above tools, between the labyrinth 
phase (or a lattice with voids) and the lattice phase (es- 
pecially in the soft-core regime when the corresponding 
RDF does not differ much from each other, and both 
phases are characterised by small values of /. However, 
these phases can be easily distinguished by employing 
another simple criterion, namely, the occupation factor 
which characterizes the fraction of the space occupied by 
particles to the total space. 

Let us assume that all the patterns are formed by is- 
lands of triangular lattice with the average distance be- 
tween two nearest neighbor particles r = ri, where t\ is 
the position of the first peak of the RDF. Then the ratio 
of the area occupied by the particles and the whole sim- 
ulation area is A = (ri/a)^. As shown in Fig. 13(a), this 
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Let us define the particles with the inter-particle distance 
shorter than ^ as neighboring particles. Then 
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FIG. 13. (Color online) The occupation factor A = (ri/a)^ 
for b = 1.1 (a). Curves A, B, C, D, E, and F correspond to 
re = 1.9, 2.1, 2.3, 2.5, 2.7, and 2.9, respectively. The occu- 
pation factor, A ior b = 4 (b). Curves A, B, C, D, E, and F 
correspond to re = 1.15, 1.3, 1.45, 1.6, 1.75, and 1.9, respec- 
tively. 



ratio can be used to distinguish the labyrinth phase and 
the lattice phase, since for the lattice phase yl « 1. 

We also found that circular clusters in the case of soft- 
core interaction are very stable with respect to an in- 
crease in the density. In this case, the additional parti- 
cles cannot efficiently increase the occupation factor due 
to the increase in the local density of the core. However, 
for the stripe and labyrinth phases, the occupation fac- 
tor A increases with the density (see Fig. 13(a)). These 
phases are not stable for high density. Thus the large 
jump in A{N) at TV = 6500 in curve C in Fig. 13(a) 
shows the phase transition from the stripe phase to the 
lattice phase. 

As one can expect, in the case of hard-core interaction, 
the occupation factor increases nearly linearly versus the 
density (see Fig. 13(b)). Therefore, the polygon shape 
and island-like clusters are not stable: with increasing 
density they evolve to the lattice phase (the plateau in 
A{N) in Fig. 13(b)). 

In addition, let us introduce another useful quantity 
which characterizes the degree of "perfection" of a lattice. 



e = 



''\ 



L 27Trg{r)n(r — a)^dr 
L 2TTrg{r)ndr 



1 

a' 



/o ''5('')(^ ~ a)^dr 
\ lo rg{r)dr 



(11) 



is an average displacement of the inter-distance between 
two neighbor particles (measured in units of a), which 
is independent of the density. The function e is non- 
zero for a deformed triangular lattice and for the ideal 
triangular lattice. Thus, e quantitatively measures the 
degree of perfection of a lattice. Note that e is only used 
for lattices as an auxiliary tool to distinguish ordered 
lattices from less ordered ones (see Table I). 



V. CONCLUSIONS 

Using molecular-dynamics simulations, we analyzed 
the pattern formation and identified different phases in a 
system of particles interacting via a non-monotonic po- 
tential, with a repulsive short-range part and attractive 
long-range part. The form of the interacting potential is 
generic: it describes, depending on specific parameters, 
the interparticle interaction in a variety of physical sys- 
tems ranging from, e.g., atoms and molecules (Lennard- 
Jones potential) to colloids and neutron stars. It can 
also be used as a model of inter-vortex interaction in 
low K type-II superconductors and in recently discovered 
so-called "type-1.5" superconductors. The obtained dif- 
ferent phases were summarized in a phase diagram in the 
plane "critical radius re — density n" (r^ is the critical 
radius where the interaction force changes its sign). We 
also analyzed the influence of the hardness of the "core" , 
i.e., the strength of the repulsive core part of the inter- 
action potential on the pattern formation. 

We developed a set of criteria in order to unambigu- 
ously identify the obtained phases using the following 
approaches: (i) the Radial Distribution Function (RDF) 
g{r), (ii) the local density function /, and (iii) the occu- 
pation factor A. In addition, we introduced a parameter 
which characterizes the degree of perfection of a lattice 
s. Employing these approaches, we elaborated the crite- 
ria for the identification of the different phases which are 
summarized in Table I. 
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TABLE I. The set of criteria used to quantitatively identify 
the different phases in terms of the RDF g{r), the local density 
function I, the occupation factor A, and the parameter e (for 
details, see the text). 
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